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The scarring effect of short unstable periodic orbits up to 
times of the order of the first recurrence is well understood. 
Much less is known, however, about what happens past this 
short-time limit. By considering the evolution of a dynam- 
ically averaged wave packet, we show that the dynamics for 
longer times is controlled by only a few related short periodic 
orbits and their interplay. 

PACS numbers: 05.45.+b, 03.65.Ge, 03.65.Sq 



The study of the quantum manifestations of classical 
chaos is at present a topic of very active research inter- 
est |lj]. Great advance came from random matrix the- 
ory (RMT) which provides an understanding of univer- 
sal statistical properties of quantum spectra |Q|. The 
most striking departure from RMT described so far in 
the literature is the phenomenon known as "scar" Q. 
This term describes an anomalous localization of quan- 
tum probability density along unstable periodic orbits 
(PO) in classically chaotic systems. Heller showed [Q the 
importance of the first recurrences of POs in the scarring 
effect. In a subsequent work Tomsovic and Heller ||] 
demonstrated that the semiclassical propagation can be 
carried out with remarkable precission long after classi- 
cal fine structure had developed on a scale much smaller 
than h, by computing the corresponding correlation func- 
tion, Csci(i), as a sum of contributions of the homoclinic 
excursions of the PO. This procedure is however cum- 
bersome, since in general many orbits are needed to ob- 
tain converged results, and this number increases rapidly 
with time (see Ref. B for details). This picture is greatly 
simplified if alternatively considering the corresponding 
averaged dynamics for finite periods of time. In this case, 
as will be shown in this Letter, a structure of a few 
short POs emerges that govern the quantum dynamics 
for times past the first recurrence of the original PO. 

Understanding scarring can be tackled from two sides. 
One way is trying to disentangle the complexity involved 
in the distribution of individual levels in the spectra of 
classically chaotic systems [^|-|[. For example, in Ref. 
[H structures localized on short POs of the stadium were 
obtained, by considering state correlation diagrams and 
iteratively removing the parametric interactions (avoided 
crossings) between the involved eigenstates. The other 
way consists of approaching the problem in a much more 



straighforward fashion, by studying how the dynamics of 
POs induce scars in the eigenfunctions of the system. 

Heller's work provided a time-dependent view that 
shows how recurrences in the short time dynamics of a 
wave packet along the neighborhood of an isolated PO 
produces the accumulation of quantum probability den- 
sity characteristic of scars. Very recently Kaplan and 
Heller [^ have shown how the use of coherent wave packet 
sums, decaying as the log-time, leads to enhanced scar- 
ring. Later, it was described by some of us ||l^ how 
(nonstationary) wave functions highly localized over POs 
can be constructed from finite time Fourier transform of 
wave packets. This method allows the calculation of the 
contribution of each eigenstate to scars. 

The time-dependent approach presents the additional 
advantage of being easily connected with spectroscopy 
experiments through Fourier transform into the energy 
domain of the autocorrelation function generated by a 
test initial wave packet [0|. In the short times limit, 
recurrences originated by isolated unstable POs in the 
vicinity of the packet translate into low resolution fea- 
tures of the corresponding spectra (or local density of 
states) consisting of peaks or bands whose widths are 
determined by their Lyapunov exponents. In the other 
limit, corresponding to long times, non-Hnear effects in- 
troduce modulations in these envelopes leading essen- 
tially to Porter-Thomas distributions of spectral inten- 
sities, although it has been shown recently ||l^ that the 
scarring effect is still noticeable in the tail of the distri- 
bution. Finally, when the Heisenberg time is reached, in- 
dividual eigenstates are resolved, and we have essentially 
the infinite resolution version of the spectrum. In this 
sense, the overall aspect of chaotic spectra is highly con- 
ditioned by the least unstable PO (LUPO) contained in 
the phase space region spanned by the test wave packet. 

In this Letter we consider the case of intermediate 
times which allow the participation of orbits longer than 
the LUPO. The associated low resolution spectra split 
into several components showing evidence of the partici- 
pation of only a few localized structures corresponding to 
short POs, which also present their signatures in the cor- 
relation function. It is the magnitude of the interaction 
(in the sense of Ref. M ) among these localized wave func- 
tions and that associated to the LUPO that determines 
which orbits come into play. 

In our study we used the desymmetrized stadium bil- 



Hard of radius r = 1 and enclosed area of 1 + 7r/4. This 
system constitutes a paradigm in the study of hard chaos. 
From the classical point of view it has been demonstrated 
to be ergodic. 

To study the scarring effect of the different stadium 
POs we closely follow the method described in Ref. [ p^ , 
which consists in following for finite times the dynam- 
ics of a wave packet initially located in the vicinity of a 
particular PO. In our case we use a harmonic oscillator 
coherent state (throughout the paper h is set equal to 
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with a=30.68. The time evolution of this state can be 
followed by projection onto the stadium eigenstates, \n): 
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which are calculated by the scaling method [|T3]. Recur- 
rences in the corresponding correlation function, C{t) = 
((/)(0)|(?!)(i)), determine the structure in the associated fi- 
nite resolution spectrum, 
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where a window function, Writ), has been introduced to 
eliminate features in C{t) taking place after a given time. 
Choosing a Gaussian form, e"* "-^ , for this window the 
following expression for It{E) is obtained 
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As shown in Ref. flQ], it is also possible to compute the 
wave functions associated to each low resolution band (or 
envelope), by means of the expression 
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These band wave functions have been shown to be highly 
localized along the scarring PO, when the time of the first 
recurrence is used as the smoothing time |1^] . 

We will concentrate our study in one of the most 
representative short PO of the desymmetrized stadium; 
namely, that running along the diagonal joining the two 
extreme points on both axis, which corresponds to a dia- 
mond shape in the full version of the stadium. This PO 
is presented in the inset to Fig. |l| along with some others 
which are relevant to our study. 
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FIG. 1. Autocorrelation function corresponding to a wave 
packet initially centered on the diagonal periodic orbit (a) of a 
desymmetrized stadium billiard with r = 1 and area 1 + 7r/4. 
This and other orbits relevant to our study are shown in the 
inset. Arrows indicate the recurrence times of these periodic 
orbits. 

Using a wave packet [Eq. (n^)] initially located on 
the middle of the diagonal PO [(xq, yo, -P?, ^°) = 
(1, l/2,96/-\/5, — 48/-y5)], the autocorrelation function 
and the corresponding infinite (T = oo) and low resolu- 
tion (T = 0.04) spectra have been computed. The results 
are presented in Figs. |l| and ^, respectively. Notice that 
the time used in the smoothed version of the spectrum 
is roughly equal to the period of appearance of the first 
recurrence in C(i), to which only the diagonal PO con- 
tributes. It can be observed that the peaks correspond- 
ing to /oo appear grouped in clumps, so that It exhibit 
a series of (seventeen) very well defined bands, which are 
equally spaced. Moreover, the positions of these peaks, 
calculated as the mean energy value of each band, are 
in very good agreement with the energies obtained from 
the wave numbers, fc, quantized according to the usual 
Bohr-Sommerfeld quantization rule: 
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where n is the number of nodes along the orbit, L = 2-\/5 
its lenght, and v = % the corresponding Maslov index. As 
an illustration, we show in the upper part of Fig. the 
wave functions (squared) associated to bands 9 to 12, 
calculated using Eq. (ph. The values of \/{E) are also 
given below each plot, so that the agreement with the 
value predicted by the quantization condition (Bh can be 
checked. As expected, the wave functions appear highly 
localized along the diagonal PO, corresponding to a num- 
ber of nodes of 31 to 34, respectively. This result could 
have also been obtained with the more elaborated coher- 
ent wave packets sums of Kaplan and Heller H . 
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FIG. 2. Infinite resolution spectrum (sticks) and smoothed 
version of it for T — 0.04 (full line) corresponding to the 
autocorrelation function of Fig. m In the upper part, the 
wave functions associated to bands 9 to 12, calculated using 
Eq. (H) and T — 0.04, are presented. The corresponding 
values of wJE) are shown below each plot. 

Let us consider now recurrences taking place at longer 
times in the correlation function of Fig. pi. One can forsee 
that their analysis will be more complicated; for one 
thing they appear at longer times so that the packet has 
had the opportunity to explore more extended regions of 
phase space (where the linearized dynamics around the 
PO is no longer valid), and also several orbits contribute 
to them. The mechanism responsible for the dynamical 
coupling of these POs, accessed by the phase space sam- 
pling associated to the choice of the initial wave packet 
(H), is not obvious. Some preliminary results [H indicate 
that the dynamics along the manifolds emanating from 
them and their proximity, is an important factor in this 
issue. 

For this purpose, let us examine the group of peaks in 
C{t) at t ~ 0.08 — 0.12. It presents a complicated struc- 
ture; in addition to the main maximum at t ~ 0.09, it 
exhibits a shoulder at smaller time values and four much 
shorter peaks at longer times. Taking into account the 
periods of the orbits presented in Fig. n^ a tentative as- 
signment of the contribution to each peak has been made. 
The result is indicated by labeled arrows over the curve 
of C{t). To check our hypothesis, we consider a series 
of low resolution spectra computed using progressively 
larger smoothing times, so that POs with increasingly 
longer periods are allowed to come into play. Some rep- 
resentative results, corresponding to the 9-th band of 
the spectrum of Fig. g, are presented in Fig. g|. Our cal- 
culations show that in the range T = 0.02 — 0.07 the 
smoothed band consists of only one peak, which begins 
to develop an incipient intraband structure as T increases 
(see different curves in Fig. g) . After that time the struc- 
ture gets more clearly defined, so that for the interval 
T = 0.12 — 0.20 four intraband components are readily 
observable (dotted line in Fig. |^). The associated wave 
functions are shown in the upper part of the figure. 
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FIG. 3. Enlargement of the region of the spectrum of Fig. 
corresponding to the 9-th band. The smoothed version has 
been computed for three different values of the resolution 
time, T: 0.02 (fuU line), 0.04 (dashed line), 0.07 (dot-dashed 
line) and 0.135 (dotted line). In the upper part, the wave 
functions associated to the labeled bands are presented. The 
corresponding values of ^ {E) are shown below each plot. 

As can be seen, the wave function labeled by 2, cor- 
responding to T = 0.07, is very similar to that already 
presented in Fig. g that was computed at T = 0.04. How- 
ever, function 1, calculated at a smaller time, presents a 
much stronger localization along the diagonal PO, being 
better defined specially at the corners. This is due to 
the fact that a smaller resolution time allows the partic- 
ipation of more eigenstates in the expansion of the wave 
function associated to each band. Finally, at the largest 
time considered, T — 0.135, the band splits into four 
components, whose associated wave functions are shown 
in plots 3 to 6 of Fig. ||. By visual inspection, these func- 
tions are found to be scarred by POs (d), (c), (b)/(a) and 
(d), respectively. This assignment was confirmed more 
quantitatively in the following way. In the first place, we 
constructed localized wave functions on the scarring or- 
bits in the same way described above for the diagonal PO 
(see Fig. ^, and calculated the overlaps of the resulting 
functions with those of the intrabands considered. The 
wave function associated to peak number 5 is a linear 
combination (56% and 42%) of the structures localized 
on orbits (b) and (a), being the corresponding overlap 
values 0.75 and 0.65, respectively. Values greater than 
0.7 were obtained for the other three peaks. In the sec- 
ond place, the corresponding calculated values of \/{E): 
45.959, 46.352, 46.809, and 47.267, agree remarkably well 
with the predicted values of k that are obtained from the 
Bohr-Sommerfeld quantization rule [Eq. (m with the ap- 
propriate values for the parameters L and v]: 45.966, 
46.367, 46.748/46.715, and 47.225. The reason why the 
9-th band in the spectrum of Fig. || splits into these com- 



ponents can be understood by considering the Hamilto- 
nian matrix elements among localized states 1 and 3-6 
of Fig. II These numbers, defined as interactions in Ref. 
|g[, when calculated (226, 664, 1706 and 533 respectively) 
turn out to be larger than any other interaction with lo- 
calized structures in this energy region. The picture that 
emerges is analogous to that of non-overlapping reso- 
nances embedded in a continuum (see discussion in Ref. 
p^). The band or intraband structures obtained with 
Eq. (||) correspond, loosely speaking, to states dynam- 
ically averaged over classical paths. They are initially 
located in the neighbourhood of unstable POs (equiv- 
alent to a decaying resonance) which are the only ren- 
mant of order in the middle of the surrounding chaotic 
sea (playing the role of the continuum). Moreover, the 
corresponding spectra present a great resemblance with 
resonance spectra, and the rhs of Eq. (||) can be viewed as 
an approximate projection operator (on each resonance) 
acting on \4>{0))- 

The same analysis can be applied to investigate the 
intraband structure of other bands in the spectrum of 
Fig. g, in order to discover the signature of other POs. 
This is the case, for example, of the wave functions pre- 
sented in Fig. H which are scarred by POs (e) and (f), 
respectively. They correspond to peaks in the intraband 
structure of bands 12 and 13 of Fig. ||. Again the com- 
puted values of \/{E): 50.282 and 52.757, agree very well 
with the predicted values of k: 50.245 and 52.722. 
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FIG. 4. Wave functions corresponding to different compo- 
nents of bands 12 and 13 in the spectra of Fig. ti, showing 
localization along periodic orbits (e) and (f) of Fig. hi The 
corresponding values of y/ (E) are shown below each plot. 

This result indicates that the intraband structure asso- 
ciated to different bands can be quite distinct. As shown, 
this is due to the fact that all POs involved in a given 
intraband pattern (corresponding to a given smoothing 
time) quantize in different ways, so that each one gives 
rise to a characteristic spacing. It is the relation among 
the values of these spacings and that of the original orbit 
(the LUPO in our case), which determines the distribu- 
tion of new peaks in each band. Since they are in prin- 
ciple unconnected, the pattern under each band will be 
different. 

In conclusion, we have presented a study of the scarring 
effect of POs beyond the short time dynamics in a very 
chaotic system. A few POs, other than the LUPO, have 
been shown to be able to develop their influence in the 
low resolution features of the corresponding spectra or 



local density of states, appearing at the associated quan- 
tized energy values; in other words they are responsible 
for the averaged dynamics in this time range. It has also 
been described how localized wave functions can be con- 
structed for this intermediate time domain. These results 
provide a good smoothed picture around a given PO and 
the family of orbits connected to it, and then constitute 
an important first step in order to fully disentangle the 
complexity involved in the eigenvalues spectrum of very 
chaotic systems. 
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